Role of the tetraquark in the chiral phase transition 
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We investigate the implications of a light tetraquark field on chiral symmetry restoration at 
nonzero temperature within a simple chirally symmetric model. In order for the chiral phase transi- 
tion to be crossover, as shown by lattice QCD studies, a strong mixing between scalar quarkonium 
and tetraquark fields is required. This leads to a light (~ 0.4 GeV), predominantly tetraquark state, 
and a heavy (~ 1.2 GeV), predominantly quarkonium state in the vacuum, in accordance with re- 
cently advocated interpretations of spectroscopy data. The mixing even increases with temperature 
and leads to an interchange of the roles of the originally heavy, predominantly quarkonium state 
and the originally light, predominantly tetraquark state. Then, as expected, the scalar quarkonium 
is a light state when becoming degenerate in mass with the pion as chiral symmetry is restored at 
nonzero temperature. 
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Introduction- In the last 30 years theoretical and ex- 
perimental work on the light scalar mesons with mass 
below ~ 1.8 GeV initiated an intense debate about their 
nature. The issue is that too many scalar resonances 
have been identified than can be accommodated in a 
naive quark-antiquark picture. For instance, in the scalar 
isoscalar channel there are five states: Jo (600), /o(980), 
/o(1370), / (1500), and / (1710) [J. In order to ex- 
plain their nature, quarkonia, tetraquark and mesonic 
molecular assignments, as well as a scalar glueball state 
with mass around 1.5 GeV suggested by lattice studies of 
quantum chromodynamics (QCD) 2], have been investi- 
gated in a variety of combinations and mixing patterns 
[EEL [Ell]- Yet, a clear answer to the question which reso- 
nance should be dominantly identified as a scalar quarko- 
nium (i.e., quark-antiquark) state is not at hand. 

Nowadays evidence for a full nonet of scalars with 
mass below 1 GeV is mounting; these are the already 
mentioned isoscalars /o(600) and /o(980), as well as the 
isotriplet a (980), and the two isodoublets of A'q(800). 
As proposed long ago by Jaffe Q, a tetraquark assign- 
ment for these states can explain some puzzling prop- 
erties, such as the mass ordering which is reversed com- 
pared to the expectation from a quark-antiquark picture, 
and the strong coupling of a (980) and /o(980) to kaons 
Q. Within this context the lightest scalar resonance 
/o(600) is interpreted as a tetraquark state g [u, d] [u, d] 
where the commutator indicates an antisymmetric flavor 
(as well as color) configuration of the diquark. Further in- 
dications of a non-quarkonium nature of the scalar states 
below 1 GeV, and thus of / (600), are obtained from a 
large- N c study in the framework of unitarized chiral per- 
turbation theory Q and in the lattice studies of Ref. [l(| • 

If the light scalars are (predominantly) tetraquark 
states, the question is how to identify the lightest quarko- 
nium state Tin = 1/ V2(uu + dd) : the broad resonance 
Jo (1370) is the first candidate. This assignment is also 
supported by the fact that the scalar quarkonia are p- 



wave states, and thus expected to lie above 1 GeV to- 
gether with other p-wave quarkonia such as axial-vector 
and tensor mesons. The two isoscalars of the quarko- 
nia nonet and the scalar glueball can mix, forming the 
states /o(1370), /o(1500), and fp(1710); such scenarios 
have been discussed in Refs. [11(. While it is not clear if 
/o(1500) or / (1710) carries the largest glueball amount, 
all the above cited works agree in the assignment of a 
dominant nn component to the resonance /o(1370). 

For Nf massless quark flavors, the QCD Lagrangian 
[including effects from the U(1)a anomaly [12|] has a chi- 
ral SU(N f ) r x SU(N f )i x U(l) v symmetry, V = r + l. If 
this symmetry is linearly realized, the mass eigenstates 
of QCD come in degenerate pairs, so-called chiral part- 
ners, which have the same quantum numbers except for 
parity and G-parity. For instance, the chiral partners of 
the scalar isoscalar meson are the pseudoscalar isotriplet 
mesons, i.e., the pions. The chiral symmetry is spon- 
taneously broken to SU(Nf)v in the vacuum [13], thus 
lifting the degeneracy of the chiral partners and render- 
ing the pions Goldstone bosons. Since the pion is com- 
monly regarded to be a quark-antiquark state, and if 
the resonance /o(1370) is predominantly a quarkonium 
state, the latter should be considered as the chiral part- 
ner of the pion, and not /o(600), as is usually assumed. 
Consequently, if the chiral symmetry of QCD is restored 
above some critical temperature T c , as predicted by lat- 
tice QCD calculations [lj], the resonance /o(1370), and 
not /o(600), should become degenerate in mass with the 
pion. Clarifying this issue is important for the interpreta- 
tion of data from heavy-ion-collision experiments, whose 
major goal is to identify signatures for chiral symmetry 
restoration at nonzero temperature T. 

The aim of this paper is a first step towards investi- 
gating this scenario of chiral symmetry restoration. We 
employ the toy model discussed in Ref. [f|, which is the 
Nf = 2 limit of a more general chiral Lagrangian for 
Nf = 3. This model contains only one tetraquark field 
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in addition to the scalar quarkonium field and the pions. 
Although the other mentioned scalar-isoscalar resonances 
/o(980), / (1500), and / (1710) are not included at the 
present stage, this model has all the essential features to 
analyze the role of the tetraquark and its mixing with 
the quarkonium at nonzero temperature. To this end we 
emp loy the Cornwall- Jackiw-Tomboulis (CJT) formalism 
[l5| | in the Hartree-Fock approximation 16]. We shall 
show that, as T increases, the /o(1370) becomes lighter 
and its tetraquark admixture grows. Within our model 
calculation, we also find that, for a large range of param- 
eters, there exists a certain temperature T s < T c above 
which the state that is predominantly quarkonium be- 
comes lighter than the state that is mostly tetraquark. At 
and above T c , the state which is predominantly quarko- 
nium becomes degenerate with the pion, as expected for 
chiral symmetry restoration. 

The model- We consider the pion triplet ff, the bare 
quarkonium field tp = nn and the bare tetraquark field 
X = \\u, d] [u, d]. The potential defining our model 
emerges as the SU(2) r x SU(2)i limit of a more gen- 
eral SU (3) r x SU (3); chiral invariant Lagrangian studied 
in Ref. [6J and reads explicitly: 

V = ^+^-Fy-s^+ 1 -My~ 9X (^ +v? 2 ) , (1) 

where e parametrizes explicit chiral symmetry breaking 
by nonzero quark masses and g, chosen to be > 0, is the 
interaction strength of the tetraquark field x [which is a 
singlet under SU (2) r x SU(2)i] with the quarkonia fields. 
As we shall see, g also determines the mixing of the scalar 
fields. When j^Oa simple linear sigma model for tp and 
7r is left. In fact, the field x, with mass M x , decouples in 
this limit. The minimum of the potential (JTJ) is, to order 
0(e), assumed for 

F e g 2 

yjl - 2.g 2 /(AA/2) 2XF 2 M* 

and 7? = 0. The nn condensate ipo is identified with the 
pion decay constant f n = 92.4 MeV. Note that the vac- 
uum expectation value (vev) xo is proportional to p\. 
Thus, the tetraquark condensate xo is induced by the 
spontaneous symmetry breaking in the quarkonium sec- 
tor. After shifting the fields tp — > tp a + tp and x ~~ * Xo + X 
and expanding the potential around the minimum, we 
obtain up to second order in the fields 

(3) 

where M 2 = p> 2 ($\ - §Q - \F 2 , A/ 2 = ^. The value 

M v = 0.139 GeV is used. Due to the off-diagonal terms 
in the mass matrix in Eq. ([3]), the fields tp and x are n °t 
mass eigenstates of the potential V. The latter, denoted 
by (H,S), are obtained after an SO(2) rotation of the 



fields (tp, x) 

( H \ _ ( cos6> sin6» \ ( x \ (£] 
\S J \ -sin0 o cos6» J \ tp J ' [ } 

where 9q — | arctan jjrr^Ti ■ The tree-level masses of H 
and S are: 

M% = M 2 cos 2 6» + M 2 sin 2 6 - 2gp Q sin(26» ), 
Ml = M 2 cos 2 6» + M 2 sin 2 O + 2gp Q sin(26» ). 

Assuming — ir/4 < 6q < 7r/4, the state H is then predom- 
inantly tetraquark and S predominantly quarkonium. As 
discussed above, we shall identify the state H with the 
resonance / (600) and the state S with / (1370). A nat- 
ural choice is then that the pure tetraquark should be 
lighter than the pure quarkonium, i.e., M x < M v . Using 
the fact that the trace and the determinant of the mass 
matrices before and after the 5*0(2) rotation are equal, 

we obtain (M| - M 2 H ) 2 = {Ml - M 2 ) 2 + {Agtpof , im- 
plying that for g > the masses of the states H and 
S repel each other: M H < M x < M v < M s and that 
\M%-M%\>4g<po. 

Note that we refrain from a more elaborate study of 
other vacuum properties of the model (such as decay 
widths), since a realistic description of the latter requires 
the inclusion of other scalar states and of (axial-)vector 
mesons [Tt|. [l8|. fl9] . 

As a side remark, we mention that the tree-level de- 
cay width of /o(600) is larger than 300 MeV when the 
mass lies above 0.6 GeV and when the mixing is large. 
However, we refrain from a more elaborate study of vac- 
uum properties, since a realistic description of the latter 
requires the inclusion of other scalar states and of (axial- 
hector mesons [l?], OH Oil and we concentrate on the 
behavior at nonzero T. 

Results- In order to study chiral symmetry restora- 
tion at nonzero T we employ the CJT-formalism in the 
Hartree-Fock approximation; for details see Ref. [l6| . As 
a result, the masses M H (T), M S (T), M W (T), and the 
mixing angle 9(T) become functions of T. Moreover, 
both scalar-isoscalar fields attain T-dependent vacuum 
expectation values (vev's), tp — > <p(T) for the quarko- 
nium and xo - > x{T) for the tetraquark state, respec- 
tively, with tp(0) = tp and x(0) = Xo, see Eq. ((2]). 

We first study the order of the chiral phase transi- 
tion and the associated T c as a function of the model 
parameters g, M H = M H (0), and M s = M s (0). Fig- 
ure 1(a) shows the phase diagram in the g-M$ plane for 
fixed Mh = 0.4 GeV which is close to the value of Ref. 
[20| . while Fig. 1(b) depicts the g-Mn plane for fixed 
Mg = 1.2 GeV which is in the experimentally established 
range of values for /o (1370) [lj. 

One observes in Fig. 1(a) that, in the limit g — * in 
which S is a pure quarkonium, the transition is crossover 
below a value Ms — 0.948 GeV and of first-order above it. 
The fact that a heavy chiral partner of the pion induces 
a first-order chiral phase transition has been discussed 
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FIG. 1: Order of the phase transition as function of the pa- 
rameters of the model. The forbidden area violates the con- 
straint \M% — Mh\ > igipo- On the border line between the 
first-order and the crossover transitions a second-order phase 
transition is realized. 
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previously, see e.g. Ref. [HI]. This means that, in a linear 
sigma model without tetraquark degrees of freedom, a 
heavy (i.e., mass larger 1 GeV) chiral partner of the pion 
is excluded by lattice QCD calculations [HI, which in- 
dicate a crossover transition. Including tetraquarks dra- 
matically changes this conclusion: as shown in Fig. 1 (a) , 
for increasing g the region of crossover transitions extends 
towards larger values of Mg. The scenario outlined in the 
Introduction, in which H = /o(600) and S = / (1370), 
can accommodate for a crossover transition if g, and thus 
the mixing between quarkonium and tetraquark, is suf- 
ficiently large. Note, however, that the range of g val- 
ues for which the transition is crossover, narrows sub- 
stantially when Ms increases. We also find that, along 
the line of second-order phase transitions in Fig. 1 (a), 
T c sizeably decreases for increasing g, for instance from 
T c ~ 241 MeV at g = to T c ~ 186 MeV at g = 3 GeV 
and T c ~ 173 MeV at g = 4 GeV. We observe in Fig. 1(b) 
that a crossover transition occurs only for small values of 
Mh- The crossover region widens when g increases. In 
order to accommodate a value ~ 0.4 GeV [2(J, a large 
value of g is required. Along the line of second-order 
transitions, T c first decreases, and then increases for in- 
creasing g. The minimum T c ~ 145 MeV occurs for g ~ 2 
GeV. 

We now study the T-dependence of masses, conden- 
sates (2lj . and the mixing angle in more detail in the 
case of M H = 0.4 GeV and M s = 1.2 GeV (in the range 
quoted by Refs. EL&HJ; a mass M H ~ 0.4 GeV, altough 
leading to a too small tree-level decay width due to lack 
of phase space, allows for a nice illustrative description of 
the qualitative features of the nonzero T). Also, we set 
the coupling strength g = 3.4 GeV, in order to obtain a 
crossover phase transition in agreement with lattice QCD 
calculations [1J|. These parameter values lead (together 
with ip = /„. = 92 A MeV and M w = 0.139 GeV) to 
M x = 0.82 GeV, M v = 0.96 GeV, and A = 52.85. 

The condensates ip(T) and x(T) are shown in Fig. 2(a). 



FIG. 2: Condensates (panel a), masses (panel b) and mixing 
angle (panel c) as function of T. 



The quark condensate f(T) drops at T c and then ap- 
proaches zero, signaling the restoration of chiral symme- 
try. The tetraquark condensate x{T) first drops together 
with <p(T), but increases above T c : this is due to the fact 
that in the equation determining x(T), the growth of the 
S and H tadpole contributions with T has to be balanced 
by an increase of x(T)\ this could be different if we in- 
clude additional terms ~ x i n the original potential |T]) . 
In any case, the increase of x(T) affects the behavior of 
the masses or of other physical quantities only slightly. 
Note that the field x is a singlet under chiral transforma- 
tions in the Nf = 2 case and therefore the nonzero value 
of the condensate does not imply a breaking of chiral 
symmetry at high T. 

The T-dependent masses M H (T), M S (T), M V (T) of 
the particles H, S, and it are shown in Fig. 2(b). The 
solid line corresponds to Mg(T), the mass of the state 
which is predominantly quarkonium (|0(T)| < 7t/4), and 
the dotted line to Mjj(T), the mass of the state which 
is predominantly tetraquark. At T s ~ 160 MeV, both 
masses behave discontinuously and the states interchange 
their roles: for T < T s , the state S is the heavier scalar 
and H the lighter one, and for T > T s , the state S is 
lighter than H. Above T s the state S becomes degenerate 
with the pions as in the sigma model without tetraquark. 
Note that, before becoming degenerate with the pion, the 
thermal mass of the lightest state (identified with H for 
T <T S and with S above it) slightly decreases, see also 
Ref. P3- 

In Fig. 2(c) the mixing angle Q(T) is plotted. At T a , 
0(T) is discontinuous, which leads to the discontinuity in 
the masses noted above: it jumps suddenly from 7r/4 to 
-7t/4, lim T _^ T ± 9(T) = =p-J. At T s the mixing is maximal: 
the two physical states H and S have the same amount 
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(50%) of quarkonium and tetraquark. Note that, at T s , 
the field tp and \ are degenerate in mass, which explains 
the maximal value of the mixing angle. As a last remark 
we note that the relative magnitude of T s and T c (T s < T c 
as in our example or vice versa) depends on the choice 
of the parameters. When increasing the mixing strength 
g the temperature T s decreases faster than T c , thus real- 
izing the ordering T s < T c in which the jump occurs at 
smaller temperatures than chiral symmetry restoration. 

Conclusions- In this Letter we proposed a novel sce- 
nario for chiral symmetry restoration at nonzero T, in 
which two scalar-isoscalar states, a tetraquark and bare 
quarkonium field, are considered. The mixing of the lat- 
ter two generates two physical states which can be associ- 
ated with the resonances /o(600) and /o(1370). When the 
tetraquark mass is smaller than the quarkonium mass, 
as supported by various spectroscopic studies in the vac- 
uum, the state Jo (600) is predominantly tetraquark and 
/o(1370) is predominantly quarkonium. This scenario 
has been studied by employing a simple model which 
includes only these two scalar resonances and the pion 
triplet. 

A remarkable aspect of our results is that the 
tetraquark-quarkonium mixing generates a softer first- 
order phase transition or, depending on the coupling 
strength, even a crossover transition. While in the stan- 
dard linear sigma model (g = 0) a heavy chiral partner 
of the pion (with mass exceeding 1 GeV) always leads to 
a first-order phase transition, this is not necessarily the 
case when tetraquark-quarkonium mixing is considered: 



for sufficiently large coupling strength g, the chiral transi- 
tion is crossover, just like in lattice QCD studies [l4| . We 
also demonstrated that the mixing between quarkonium 
and tetraquark states increases with increasing temper- 
ature, and, in most cases, reaches its maximal value of 
45° at a temperature T s where the physical states consist 
of an equal amount of quarkonium and tetraquark. For 
T > T s the physical states interchange their roles: the 
lighter state is predominantly quarkonium and the heav- 
ier predominantly tetraquark. Further increasing T leads 
to the standard scenario of chiral symmetry restoration, 
where the scalar quarkonium becomes degenerate in mass 
with its chiral partner, the pion. Thus, our approach can 
possibly solve an inconsistency between low-energy spec- 
troscopy, where a non-quarkonium structure for /o(600) 
is favoured, and studies at nonzero temperature, where 
the scalar partner of the pion should be sufficiently light 
(~ 0.6 GeV) in order for the chiral symmetry restoring 
transition to be crossover. 

Since the present work is a first explorative study on 
the relevance of the tetraquark at nonzero T, we omit- 
ted other scalar-isoscalar states such as /o(980), /o(1500), 
and /q(1710), which would naturally appear in an SU (3)- 
symmetric model with two scalar nonets and a glueball 
state. Also (axial-)vector mesons should be considered 
[l7l IHi ] . All these fields are important in a more realis- 
tic framework which aims to describe at the same time 
vacuum phenomenology and nonzero temperature prop- 
erties. We regard the results of this paper as a motivation 
to undertake this more ambitious step in the near future. 
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